Title: Flexible Estimation of Groundwater Service Values and Time Preferences
Authors: Grant West, Heather Snell, Kent Kovacs, Rodolfo Nayga Jr.
Date: Jan 4, 2021

For questions regarding the coding and files, please e-mail Heather Snell (haprice@uark.edu).

Description: The three folders contain the files necessary to reproduce all
tables in the aforementioned article for JAERE.  MATLAB version R2018b was used
for model estimation.

> DLML
  > General
    'FlexibleRateBestSupport.m' is the general script used for creating Table 4
    and the tables in appendix B. 'FlexibleRateThumb.m' is the script used for creating
    appendix table A2.

    In either matlab script, the global variable, 'discountfun', can be set to
    "QuasiHyper", "Exponential", "Harvey", or "HM".  The diary name is currently
    set to "QH.out", but could be modified. The diary creates an "out" file which is
    essentially a text file with the output of the console.  The 'filename' determines
    what name the workspace is given when it is saved, and the
    'modelname' determines what name is used for the figure png file.

    To run bootstrapping, 'WantHessian' should be changed to 0 and 'WantBoot' should
    be set changed to 1.

    All other matlab code files support the estimation of the flexible mixing
    models with discounting and should not be changed.

    The matlab data file is 'discouting.data.mat', and the csv files (infra.time,
    jobs.time, price.time, quality.time, SW.time, wlife.time, CAP.time, and
    buffer.time) contain the values for the attributes for each time period and
    each choice question. These should not be changed.


  > ANA
    'FlexibleRateANA.m' is the script for generating table 7 (model II.b ANA).
    Other matlab files support the estimation.  These files are currently only
    set to run the quasi-hyperbolic discounting function.

  > Consequentiality
    'FlexibleRateANA.m' is the script for generating table 7 (model II.a Consequentiality).
    Other matlab files support the estimation.  These files are currently only
    set to run the quasi-hyperbolic discounting function.

> DMIXL
  'dmixl.m' is the general script used for creating appendix tables A1.  It runs the
  mixed multinomial logit models with discounting.

  In this script, the global variable, 'discountfun' is currently set to "Exponential"
  but can be set to "QuasiHyper", "Harvey", or "HM".  The diary name is currently
  set to "dmixl.exp.out", but could be modified. The diary creates an "out" file which is
  essentially a text file with the output of the console.  The files 'dmixl.exponential.out',
  'dmixl.harvey.out', 'dmixl.hm.out', and 'dmixl.qh.out' are the results used in the appendix.
  The 'modelname' determines what name the workspace is given when it is saved and
  is currently set to 'dmixl.exp.mat'.

  To reproduce the results in table A1, the appropriate start values for B and W
  (lines 149 and 150) should be set as in the excel file "DMIXL_Start_Values.xlsx"
  according to the discounting type.

  All other matlab code files in this folder support the estimation of the mixed
  models with discounting and should not be changed.

  The matlab data file is 'discouting.data.mat', and the csv files (infra.time,
  jobs.time, price.time, quality.time, SW.time, wlife.time, CAP.time, and
  buffer.time) contain the values for the attributes for each time period and
  each choice question. These should not be changed.
